Quantum-correlated motion and heralded entanglement of distant optomechanically 

coupled objects 
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The motion of two distant trapped particles or mechanical oscillators can be strongly coupled by 
light modes in a high finesse optical resonator. In a two mode ring cavity geometry, trapping, cool- 
ing and coupling is implemented by the same modes. While the cosine mode provides for trapping, 
the sine mode facilitates ground state cooling and mediates non-local interactions. For classical 
point particles the centre-of-mass mode is strongly damped and the individual momenta get anti- 
correlated. Surprisingly, quantum fluctuations induce the opposite effect of positively-correlated 
particle motion, which close to zero temperature generates entanglement. The non-classical correla- 
tions and entanglement are dissipation-induced and particularly strong after detection of a scattered 
photon in the sine mode. This allows for heralded entanglement by post-selection. Entanglement 
is concurrent with squeezing of the particle distance and relative momenta while the centre-of-mass 
observables acquires larger uncertainties. 

PACS numbers: 37.30.+i, 03.65.Ud 



I. INTRODUCTION 

The past decades have seen tremendous success in the 
implementation of control schemes for the motional state 
of matter via light fields either in free space or in optical 
cavities. A diversity of examples exist where the quan- 
tum regime of motion has been reached. The masses span 
many orders of magnitude, from the microscopic atomic 
size systems such as atoms in optical cavities [1-4] and 
laser-cooled ions in ion traps [5] to the macroscopic level 
with cavity-embedded membranes [6], mirrors [7] or lev- 
itated dielectric nano-particles [8]. 

A common interaction Hamiltonian that well ap- 
proximates many quantum light-matter interfaces is 
quadrature-quadrature coupling [9]; more specifically, 
the displacement of the mechanics is coupled directly to 
a quadrature of the high-Q optical field mode that can 
be then used as an observable for indirect position detec- 
tion. Adding a second mechanical system coupled to the 
field then allows one to engineer an effective two-particle 
mechanical coupling by eliminating the mediating light 
mode. Recently, an expansion to quadratic coupling has 
been proposed [10] and the investigation of dissipation- 
induced [11-13], noise-induced [14] and remote entan- 
glement [15, 16] has been of great interest, including a 
scheme for sensitive force measurements [17] and entan- 
glement of macroscopic oscillators [18, 19]. 

Here we show that all this can be implemented in a 
system consisting of two particles strongly trapped in the 
cosine mode of a ring cavity, where the two-particle inter- 
action is carried by sideband photons in the sine mode. 



For deep trapping it yields the typical linearized optome- 
chanical Hamiltonian [20]. First we present the general 
model of two particles moving within a ring resonator. 
We then analyse single quantum trajectories depicting 
strong correlations and entanglement. A subsequent in- 
vestigation of momentum correlations reveals classically 
forbidden positive values in steady state, even in the ab- 
sence of entanglement. The steady state shows a strong 
derealization of the centre-of-mass independent of the 
particle separation. We also show how to generate en- 
tanglement either in a pulsed regime or heralded by the 
detection of photons. Analytical calculations are car- 
ried out in the regime of strong particle confinement and 
matched to the more generally valid numerical simula- 
tions with good agreement. Finally, the occurrence of 
correlations in the system is explained in a simple adia- 
batic model. 
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We structure this presentation in two major parts ac- 
cording to the Gaussian/non-Gaussian nature of the sys- 
tem dynamics. Entanglement on trajectories and that 
on average as well as momentum correlations during the 
no n- Gaussian dynamics of the two atoms in the cavity 
are numerically investigated in the beginning, while the 
Gaussian approximation makes the topic of the second 
major part of the paper. This part corresponds either to 
atoms in very deep traps (to compare to the results ob- 
tained in the first part) or more generally to a wide range 
of optomechanical systems in the linear regime and with 
diverse physical realizations (moving membranes, levi- 
tated nano-spheres etc.). 
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FIG. 1: Sketch of the system. The pumped standing- wave 
mode (red) traps the particles, the second orthogonal one 
(black) mediates an effective interaction between them. 



II. MOTION OF TWO PARTICLES IN A RING 
RESONATOR 

We study two small polarizable particles confined 
within a symmetrically-pumped ring resonator, see fig- 
ure 1 . Symmetric pumping results in a standing- wave op- 
tical potential with spatial dependence cos 2 (kx) [20, 21]. 
The cosine mode is strongly pumped and approximated 
by a highly excited coherent state \a c ) with \a c \ ^> 1 (and 
without loss of generality a c G R). The particles scat- 
ter photons into the unpumped orthogonal sine mode. 
This setup can be generally described by the Hamilto- 
nian [20, 21] 



1 = 1 



— h hU^a 2 cos 2 (kxi) + HUqo) a sin 2 (kxi) 



HU a c 



(a + a 1 ") ^ sin(2fcxi) - HA c a ] a. (1) 



1=1 



Here, a denotes the annihilation operator of the 
quantum- mechanical sine mode, Uq < the optical po- 
tential depth per photon, X{ and pi the particles' centre- 
of-mass position and momentum operators, respectively, 
and m the particle mass. The pump is detuned by 
A c := co> p — uj c from the bare cavity resonance frequency 
uj c . To avoid instabilities we restrict ourselves to red 
detuned lasers (A c < 0) for which a cooling regime ex- 
ists [22]. The sine mode is only weakly populated by scat- 
tering such that (a) a) is negligible compared to a 2 [20]. 
Damping of the cavity mode is taken into account by 
the Liouvillian Cp = k {2a pa) — a) ap — pa) a) [23] in the 
master equation 



P=- H [H, P] 



Cp. 



(2) 



To get some insight into the dynamics, we solve equa- 
tion (2) numerically for some typical parameters. The 
direct solution of the master equation (2) is computa- 
tionally very demanding owing to the large Hilbert space 
of the joint particles-field system. Therefore, we re- 
sort to Monte Carlo wavefunction simulations [24], in 



which the system is coherently evolved between the so- 
called "quantum jumps". These jumps correspond to a 
photon detected at the resonator output [24]. Besides 
the more favourable usage of computer resources, single 
trajectories also provide an additional insight into the 
microscopic processes in the system. The simulations 
were efficiently implemented with the freely available [39] 
C++QED framework [25, 26] and performed in a joint 
momentum- and Fock basis. 

A typical trajectory is shown in figure 2, where the blue 
arrows indicate the times at which jumps occur. Initially, 
the particles were prepared in the ground state of two 
separated potential wells and owing to the deep poten- 
tial tunnelling is strongly suppressed. The momentum 
correlation coefficient C p for the two particles is defined 
as 



Cov(pi,p 2 ) (P1P2> - (Pl> (P2> 



Ap 1 Ap 2 



Ap 1 Ap 2 



(3) 



where C p = 1 means perfect correlation and C p = — 1 
perfect anti-correlation of the motion. Quantum jumps 
by photodetection trigger strong correlations and entan- 
glement between the particles due to the cavity-mediated 
interaction. The logarithmic negativity [27] already af- 
ter the first jump approaches the value for a maximally 
entangled Bell state and the correlation reaches a value 
of C p ~ 0.5. The emerging state corresponds to a su- 
perposition of two particles moving to the right and two 
particles moving to the left such that the centre-of-mass 
momentum remains zero. As we will see later this be- 
haviour is caused by the excitation of a single particle 
into the first excited state within the trap. Beginning 
with the second jump the system is subject to fast oscil- 
lations which are of the order of the trap frequency. Note 
that only the field is dissipative in our model and induces 
quantum jumps. However, the particles respond to the 
sudden changes of the field in a correlated way. 

Interestingly, while entanglement is quite pronounced 
on single trajectories, it remains small on average. The 
averaged steady-state logarithmic negativity is Sj\f ~ 
10 -3 for our parameters. The logarithmic negativity is 
not an expectation value and hence its value for an en- 
semble of trajectories cannot be obtained from its val- 
ues on single realizations — it has to be directly com- 
puted from the density operator. It is therefore possible 
that the logarithmic negativity of a mixture of entan- 
gled states is smaller than the weighted average of the 
individual values. A prominent example for this is the 
equal statistical mixture of two Bell states, which is not 
entangled. 

We now examine the momentum correlations in more 
detail. Classical simulations reveal much stronger damp- 
ing of the centre-of-mass motion than that of the relative 
motion. Hence the particles become anti-correlated [1, 
22, 28], see figure 3. Surprisingly, the quantum simu- 
lations of the ring resonator system yield the opposite 
result. On average, initially uncorrelated particles be- 
come positively correlated due to the cavity input noise. 
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FIG. 2: Photon number, momentum correlation and entan- 
glement for a single trajectory. The dashed vertical lines in- 
dicate quantum jumps. Parameters: a c — 150, Uo — — ^cjr, 
A c = —k and k — IOcjr. 
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FIG. 4: Diagonal elements of the two-particle reduced den- 
sity matrix in momentum space, initially (left) and in steady 
state (right). Ensemble average of 5000 trajectories. Same 
parameters as in figure 2. 
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FIG. 3: Classical simulations [22] for 1000 initial conditions. 
Left: initial condition, right: distribution at cjr£ = 3000. The 
particles become anti-correlated and cooled. Same parame- 
ters as in figure 2. 



This effect is visible in the steady-state density ma- 
trix presented in figure 4. There, the momentum dis- 
tribution is elongated into the first quadrant, which is 
a signature of positive correlations. Note that quan- 
tum mechanics allows pure states with positive correla- 
tions, but still zero average centre-of-mass motion with 
large uncertainty. The time evolution of the momen- 
tum correlation coefficient is depicted in figure 5. We 
observe smaller c orrelations a round the cooling sideband 
A c = —uj := — ^yAH\Uo\a^E^/h as compared to the other 
chosen detuning A c = —k. Here = fkjR := h 2 k 2 /2m 
is the recoil energy. For comparison, we also present 
the results of simulations containing the quadratic op- 
tomechanical Hamiltonian (6) introduced later on in this 
paper. The latter yields accurate results provided that 
the position spread remains small compared to a wave- 
length, kAx <C 1. This is fulfilled for operation near the 
cooling sideband, but not for A c = —k. Positive cor- 
relations of the particle motion can also be observed in 



FIG. 5: Red solid lines: momentum correlation coefficient ob- 
tained from the master equation containing (1) for A c = — k 
(upper curve) and A c = — uo (lower curve). Blue dashed lines: 
result of the oscillator approximation (6) (600 trajectories). 
The other parameters are the same as in figure 2. 



much shallower potentials where the particles are barely 
trapped and both light modes are treated quantum- 
mechanically [29]. 

In the dynamics, correlations and entanglement are 
generated by quantum jumps induced by photon count 
events. Hence they are particularly strong immediately 
after a jump. In real experiments, it is generally not pos- 
sible to exactly keep track of each emitted photon due 
to the detector efficiency, i.e. one cannot exactly follow 
a certain quantum trajectory. Hence, the system always 
evolves into a mixed state and if a photon is observed 
(measured) outside of the resonator, it is impossible to 
determine whether it was the first, the second and so on. 
Naturally, the question arises whether effects observed 
on single trajectories (like entanglement) "survive" this 
averaging and can still be expected to be observable after 
jumps. Every time a photon is detected, the mixed state 
inside the resonator is projected into the state 



ft 



apa } 



Tr (apaJ) 



(4) 



where p is the density matrix evolved according to the 
master equation (2). Interestingly, the momentum cor- 
relation coefficient (3) is nearly constant (C p « 0.5 for 
A c = — k and C p w 0.3 on the cooling sideband, respec- 
tively), regardless of the time the jump occurred at, see 
figure 6 also. The heralded entanglement measured by 
the logarithmic negativity is smaller than on single tra- 
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FIG. 6: Diagonal elements of the reduced two-particle condi- 
tional density matrix (4) after a photon detection in steady 
state. Particle correlations as well as entanglement become 
much more pronounced compared to the density matrix shown 
in figure 4. Same parameters and colour code as there. 



jectories, but still prominent with a value Ej^f ~ 0.2 for 
our parameters and a jump in steady state. 

The conditional density matrix (4) is closely related to 
the field autocorrelation function [23] 



(ototoo) _ < ata ) f 



{of a) 



(5) 



It has the very intuitive interpretation as the ratio be- 
tween the photon number after and prior to a jump 
in steady state [30]. For our parameters, it indicates 
photon bunching close to a thermal (chaotic) state (for 
which g( 2 \0) = 2). This is consistent with the picture of 
the mode being incoherently populated through particle 
noise — perfectly localized particles do not scatter. Hence 
the field expectation vanishes and only its variance gives 
a contribution. 



III. GAUSSIAN OPTOMECHANICAL 
TREATMENT 

So far, we have studied the mainly numerically accessi- 
ble general system (1) of two particles in a ring resonator. 
Now, we investigate the special case of tightly confined 
particles allowing for more analytical insight into the dy- 
namics. In this limit, it is justified to keep only the fir st- 
and second-order terms in the expansion of the trigono- 
metric factors and the Hamiltonian (1) can be mapped 
onto the linearized optomechanical model 



H = Y^ tiublbi - HA^a + hg ^(6* + b\)(a + a f ) (6) 



as shown in [20]. Here we have defined uo := 
<\/4:h\Uo\a%EYi/h, g := Uoa c k£o/V2, the oscillator length 

£o •= yh/muj and bi := (#i/£o + ^oPi /K) / V%- The scal- 
ing of the oscillator length suggests that the Lamb-Dicke 
regime fc£o <C 1 may also be reached for very heavy parti- 
cles in shallower traps. As we do not consider any direct 
particle-particle interactions, the Hamiltonian (6) is also 



valid for particles not confined in the same, but rather in 
distant sites within the resonator. 

Interestingly, the Hamiltonian (6) applies to a whole 
class of systems. A few well- studied realizations are (i) a 
cavity with vibrating end mirrors where the mirror-light 
interaction is always linear and governed mainly by the 
cavity length and the zero point motion of the mechan- 
ics [31, 32], (ii) two light membranes inside a cavity field 
positioned at the maximum slope of the field amplitude 
where the coupling depends on the reflective properties of 
the membranes and is increased with decreasing mass [33] 
and (hi) very light nano-sized dielectric spheres held in- 
side the cavity field either by an external trapping light 
mode or by means of optical tweezers [8]. 

The Hamiltonian (6) is quadratic in the bosonic oper- 
ators so that an initially Gaussian state remains Gaus- 
sian throughout the time evolution. Gaussian states are 
completely defined by a displacement vector and a co- 
variance matrix Vij := \ [Cov(&,£j) + Cov(£j,£;)], with 
£ := (x 1 ,p 1 ,x 2 ,p2iX,P), Xi := £;/£ = (bi + b\)/y/2, 
Pi : = Pi£o/h = (bi — b\)/(iy/2) and the field quadratures 
X and P [34]. As shown in [35], the time evolution of V 
is determined by the equation (its steady-state version is 
called the Lyapunov equation) 



V(t) = AV(t) + V(t)A T + B, 



(7) 



where A and B are the drift- and diffusion matrices 
appearing in stochastic Heisenberg-Langevin equations 
equivalent to the master equation. The steady-state so- 
lution of (7) for the covariances is Cov(pi,p2) = Var(p^) — 
1/2 and Qov(x\,X2) = Var(x^) — 1/2 (both particles 
behave the same way). It only exists in the cooling 
regime A c < since the momentum variance is only 
positive for red detuning. The momentum covariance 
is also genuinely positive. For C p we find the simple (g- 
independent) expression 



_ /. 2 + (A c +^) 2 



(8) 



which is precisely the ratio of the Stokes- (T_, heating) 
and anti-Stokes (T + , cooling) scattering rates found when 
adiabatically eliminating the cavity field [9] . These rates 

2 

r± = / ^2 + ^ / ^ u; )2 also define the time scale r ~ (T+ — 
T_) _1 on which the steady state is reached. 

C p has a minimum at A c = — V^ 2 + oj 2 and approaches 
unity for |A C | — >• oo and A c — >• 0~. From the form of 
the steady-state covariance matrix we can deduce that 
steady-state entanglement between the particles can only 
occur if Cov(xi,X2) < 0, which implies squeezing of the 
position variable. This is a direct consequence of the en- 
tanglement criterion derived in [36]. See figure 7 for an 
example of the Lyapunov time evolution. For such deep 
potentials (a regime which we can numerically analyse 
by changing the effective parameters uo and g to the val- 
ues listed in figure 7) we restricted ourselves to Monte 
Carlo simulations of the master equation containing the 
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FIG. 7: Upper plot: momentum correlation obtained from 
Monte Carlo simulations in the Lamb-Dicke limit for k£o = 
0.1 (red, ensemble of 500 trajectories) and solution of equa- 
tion (7) (blue). The position correlation is likewise pro- 
nounced. Lower plots: diagonal elements of the reduced 
particle density matrix in steady state (left) and of the re- 
duced conditional density matrix (4) (right). Parameters: 
uj = 200cc;r, g = 5cjr, k = IOOcjr and A c = — 20cjr. Note 
that here the particles-field interaction is much stronger than 
in figure 2 where g « 0.2cjr. 



Hamiltonian (6). Again, the C++QED framework pro- 
vides a helpful basis for their numerical implementation 
as it also supports the simulation of coupled oscillators. 

The steady-state covariance matrix reveals strong cor- 
relations between the particles and the field quadratures, 
making it intuitively clear that the reduced conditional 
density matrix for the particles strongly differs from its 
steady-state counterpart. Each photon detection highly 
influences the joint state of particles and cavity field. We 
show both density matrices in figure 7. The logarithmic 
negativity for the conditional density matrix is Sj\f ^0.25 
and C p is found to be C p « 0.9. Due to the particles-field 
correlations in steady state the conditional density ma- 
trix (4) does not describe a Gaussian state. Indeed, it 
reveals a double-peak structure. 

Given that steady state dynamics shows little or no 
traces of the entanglement that is clearly present in in- 
dividual trajectories, we pursue now an alternative road 
that investigates the transient regime. Such a regime 
could for example be reached with short light pulses that 
can be used to generate build-up of motional correla- 
tions in the system. While we are free to explore any 
regime numerically, the problem can be analytically tack- 
led mainly in the adiabatic case where the cavity mode 
can be eliminated from the dynamics, when g <C k or 
g <C |A C ±uj\ [9]. We verify, however, that the generality 
is not lost since numerical investigations show the opti- 
mal entanglement regime indeed being around the point 
analytically treated. 

We skip the cumbersome analytical elimination proce- 




FIG. 8: Transient entanglement (lower red curve) and mo- 
mentum variance (upper blue curve) as obtained from the 
Lyapunov equation. Black dashed lines: solution of the mas- 
ter equation for the ring resonator (2000 trajectories), for 
comparison. Parameters: uj — 30cjr, g — 5cjr, k = 5cjr 
and A c = — (uj — k). 

dure (that closely follows the one used in [9]) and simply 
make use of the simple form of the end result that shows 
displacement-displacement coupling of the reduced bi- 
partite system, ~ —hT ^2j k Xj%k, with 

/ A r - uj A c + uj \ 
T = - — T_ + c T + . 9 

y tv K J 

A second result of such an elimination procedure is that 
decay of the system can occur in a correlated fashion, 
namely with an effective jump operator proportional to 
b\ +62. Now, we choose the operating conditions such 
that X1X2 reduces to a simple beam splitter interac- 
tion b\b2 + b\b\ and point out a very general conclusion 
(see [37]) that this interaction does not lead to entangle- 
ment when starting in non-correlated initial states but 
rather to a simple state swap process where in the limit of 
strong coupling (coupling larger than decoherence rates), 
a quantum state of a subsystem can be written on the 
other subsystem. Than, taking advantage of the cor- 
related decay, we show that transient entanglement is 
present in transient dynamics. This is illustrated in fig- 
ure 8 where we show the time evolution of the logarithmic 
negativity. The relatively large entanglement obtained at 
half of the interaction time (defined by the inverse of the 
effective coupling strength g) can be accessed in an exper- 
imental situation for example by the use of laser pulses 
of tailored shape and duration. 

Our effective optomechanical Hamiltonian (6) is for- 
mally similar to a model describing the microscopic 
dynamics of self-organization of particles in a cavity 
field [38]. The behaviour of single trajectories there can 
be explained with the help of an adiabatic model. We 
now strive to develop a model in the spirit of [38], ca- 
pable of describing the behaviour of single Monte Carlo 
trajectories of the master equation containing the Hamil- 
tonian (6). To this aim, we diagonalize the scattering op- 
erator Ylifii + b\) to obtain the new particle basis {|e^)} 
containing radiative (eigenvalue A^ 7^ 0, they always ap- 
pear with both signs ±1) and non-radiative (A* = 0) 
states. Ignoring terms stemming from the oscillator en- 
ergies for the moment, the radiative states radiate a field 
\\iOt) with a := —\g/(n — iA c ). The crucial assump- 
tion now is adiabaticity (cavity decay assumed to de- 
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fine the fastest time scale), i.e. we assume each particle 
state \ ei) to be correlated with its associated unperturbed 
field state to find the approximated joint particle-field 
stochastic state vector 

|V^)> = £^)K->|A,a). (10) 

3 

The coefficients Ci (t) are determined by the effective con- 
ditional Monte Carlo time evolution generated by the 
non-Hermitian Hamiltonian (abbreviating \i) := |e*)) 

F nH - ^ E w E b"> i*>c/i - ^ E i A * a i 2 km 

ij k i 

(ii) 

for the particles only. Its Hermitian part contains the 
oscillator energies expressed in the scattering basis and 
couples the constituents of the latter. 

From the form of the state vector (10) one can conclude 
what happens if a photon is detected. Applying the jump 
operator ex a on the internal state (10) at time t- } yields 
a\tjj(tj)} oc Ci(tj)Xi \ei) \\ia). Hence, all information 
about the non-radiative contributions are erased and the 
relative phases between the radiative states change sign 
(as for each there exists also an eigenvalue — A*). Since 
|A^a| 2 <C 1 the reduced density matrix for the parti- 
cles directly after a jump mainly contains the pure state 
l^p) J2j c j(t)^j \ e j) (during the time evolution the par- 
ticle state reads |^ p ) « ^ Cj(t) \ej)). Inserting the coef- 
ficients for a state describing two particles in the ground 
state, the first jump excites one of them and the photon 
number increases by a factor of 3. The resulting state 
of the particles |^ p ) w -±= (l^o^i) + l^i^o)) is entangled 
with £j\f = 1 and correlated with C p = 1/2. Beginning 
with the second jump, fast oscillations with a frequency 
determined by uo build up, which is in qualitative agree- 
ment with single trajectories of the ring resonator model 
as shown in figure 2. This way it is possible to qualita- 
tively and quantitatively reconstruct single trajectories 



of the optomechanical model (6). 



IV. CONCLUSION AND OUTLOOK 

Starting from two particles in a ring resonator, our con- 
clusions are quite general. We have proposed a tripartite 
optomechanical system with two identical oscillators that 
can be correlated via the mediation of a photon field and 
the interaction amplified by the cavity confinement. The 
Gaussian evolution of the reduced mechanical biparti- 
tion has been followed both numerically and analytically, 
while non- Gaussian effects owing to the deviation from 
the bilinear Hamiltonian have been investigated numer- 
ically. The non-intuitive and classically forbidden posi- 
tive momentum correlations are an important result of 
the paper, and are interpreted as strong derealizations 
(superpositions) of the centre-of-mass independent of the 
particle separation and thus of the extension of the effec- 
tive system; this suggests possible use of the setup to test 
the quantum-classical boundary. Furthermore, we have 
studied entanglement on single Monte Carlo trajectories 
as well as of the averaged density matrix. Entanglement 
heralded by photodetection has also been investigated. 
Transient entanglement as shown here can also be ex- 
ploited in a pulsed regime where light pulses can drive the 
mechanical bipartition into a strongly entangled state. 
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